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. ABSTRACT 

\Q • 

Aims. We aim at deriving the excitation conditions of the interstellar gas as well as the local FUV intensities in the molecular cloud surrounding 
NGC 3603 to get a coherent picture of how the gas is energized by the central stars. 
Methods. The NANTEN2-4m submillimeter antenna is used to map the [CI] 1-0, 2-1 and CO 4-3, 7-6 lines in a 2' x 2' region around the 
, young OB cluster NGC 3603 YC. These data are combined with C ls O 2-1 data, HIRES-processed IRAS 60 yum and 100 fim maps of the FIR 
■ • continuum, and Spitzer/IRAC maps. 
, Results. The NANTEN2 observations show the presence of two molecular clumps located south-east and south-west of the cluster and 
, confirm the overall structure already found by previous CS and C ls O observations. We find a slight position offset of the peak intensity of CO 
t I ■ and [CI], and the atomic carbon appears to be further extended compared to the molecular material. We used the HIRES far-infrared dust data 
C^S to derive a map of the FUV field heating the dust. We constrain the FUV field to values of x = 3 - 6 x 10 3 in units of the Draine field across 
the clouds. Approximately 0.2 to 0.3 % of the total FUV energy is re-emitted in the [CII] 158 ycn\ cooling line observed by ISO. Applying 
LTE and escape probability calculations, we derive temperatures (7mmi = 43 K, Tumi = 47 K), column densities (ATmmi - 0.9 x 10 22 cm -2 , 
Nmm2 = 2.5 x 10 22 cirr 2 ) and densities («mmi = 3 x 10 3 cirr 3 , n MM2 = 10 3 - 10 4 cirr 3 ) for the two observed molecular clumps MM1 and MM2. 
Conclusions. The cluster is strongly interacting with the ambient molecular cloud, governing its structure and physical conditions. A stability 
analysis shows the existence of gravitationally collapsing gas clumps which should lead to star formation. Embedded IR sources have already 
been observed in the outskirts of the molecular cloud and seem to support our conclusions. 

Key words. ISM: clouds - ISM: structure - ISM: molecules - Submillimeter 



1. Introduction 



Send offprint requests to: M. Rollig, 
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* Founded by merging of the Sternwarte, Radioastronomisches 
Institut and Institut fur Astrophysik und Extraterrestrische Forschung 



Understanding the feedback between star formation and the in- 
terstellar medium (ISM) is impor tant in order to understan d the 
process of star formation (see IZinnecker & Yorkd l2007, and 
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Fig. 1. 8 /mi Spitzer/IRAC band (in MJy/sr) overlayed with 
CO 4-3 NANTEN2 observations presented in this paper. The 
position of the cluster is marked by a star. The boxes mark the 
positions of the pillar heads. The crosses mark the positions 
along the cuts into MM1 and MM2. IRAC pixel size is 1.2". 
The position (0,0) corresponds to R.A.(J2000.0)=1 l h 15 m 08 s .85, 
Dec.(J2000.0)=-61°16'5070. 



references therein). The physical and chemical properties of 
stars are the heritage of their parental clouds. Stars are born 
from interstellar gas and release metal-enriched material to the 
ISM when they die. Radiation from the stars is the prime heat- 
ing source for gas and dust in nearby molecular clouds and 
may trigger velocity and density fluctuations, that stimulate 
further star formation. To understand the formation of these 
next generation stars, it is important to understand ho w the pre- 
vious generation interacts with its parental clouds dDale et al 



2005). The energy incident on the clouds in the form of stellar 
far ultraviolet (FUV: 6 eV< hv < 13.6 eV) radiation is coun- 
tered by cooling continuum radiation plus emission of atomic 
fine-structure lines ([CII] 158 fim, [CI] 370 /mi, 610 /mi, and 
[OI]63 /vm, 145 /mi) and by molecular rotational lines (CO, 
H2O, OH, etc.). The cooling emission carries the imprint of 
the local physical and chemical conditions and can be used to 
infer the type of environment conducive to maintaining star 
formation. On their way into a molecular cloud FUV pho- 
tons are absorbed, and a depth-dependant chemical balance is 
established. Models of these so- called photon dominated re- 
gions (PDRs; see references in iHollenbach & Tielens 1999; 
Rolli g et al] 120071) typically predict a stratification of species 
like H/H2, C + /C/CO, and others. This strat ification is found to 
be ind ependent of the shape of the PDR dGorti & Hollenbachl 
20021) . However, observ ations often reveal [CII emission coin- 



cident with that of CO (Tatematsu et al. 1999; Ikeda et al 



Mookeriea et aL 20061 : ISakai et all 120061 : iKramer et alJ 



2002 



2008 



Sun et all 20081) . Plane-parallel models cannot explain this be 



havior unless full face-on orientation is assumed. Spherical 
models under isotropic illumination show well-correlated [CI] 
and CO emission for low density clouds (n < 10 4 cm ~ 3 ) but 



predi ct a limb brightening of the [CI] emission dRollig et al 
2006). In this work we analyze the distribution of atomic 



carbon and warm CO in the NGC 3603 star forming region 
and use LTE approximations a nd an escape probability model 
(Stut zki & Winnewisserj[l985l) to derive the excitation condi- 
tions of the gas as well as gas abundances. We refrain from ap- 
plying more detailed models like the KOSMA-t PDR mode l 
dStorzer et al.1 1 19961: iRollig et alJ 120061: ICubick et al.l 120081) . 
Recent Herschel/HIFI observations of massive star forming re- 
gions provide additional strict constraints on the PDR modeling 
of these regions, resulting in significantly alter ed predictions of 



the local physical and structural conditions (Ossenkopf et al 



120101: iDedes et alJl201dl) . Hence, we postpone all attempts of 
PDR modeling in NGC 3603 until the region has been observed 
with Herschel in the framework of the Herschel Guaranteed 
Time Key Project Warm And Dense ISM. 



2. The NGC 3603 star forming region 

NGC 3603 is located in the Carina spiral arm (/ * 291.6°, b ~ 
-0.5°) at a distance of ap proximately 7-8 kpc (see discus- 
sion in Melen a et al. 2008). It is one of the most luminous 



(Lboi > 10 7 L Q ), optically visible HII regions in the Galaxy, 
with the massive O B cluster NGC 3603 YC (You ng Cluster) 
as a power source ( Goss & Radhakri shnanl T%9l) . The com- 
pact core of the cluster was designated HD 97950 due to 
its star-like appearance. As a comparison, NGC 3603 is 100 
times more luminous than the Trapezium cluster in Orion. 
This young cluster can be considered one of the few Galactic 
starburst clusters, which are essential for the understanding 
of extragalactic star bursts. With a dynamic s tellar cluster 



mass of 17600+3800 M (Ro chauetalJ 120101) residing in 



a clo ud of total gas mass of 4 x 10 5 M G dGrabelskv et al 



19881) . it is the most compact Galactic star-forming com 



plex outside the Galactic center region (Stol te et al.l 12004). 
In a number of recent studies, the stellar cluster and 
the sur rounding HII reg i on has been investigated thor- 



eg io 

oughlv dPandev et alJl2000HSung & Bessel]||2004t lStolte et al 



Haravama et alJ 2008 



2006; Lebouteiller et al. 

iNiirnbergeri 120081 ICrowther et al.ll2010l) . To the south of the 
cluster we find a giant molecular cloud (see Figure Q}. The 
intensive radiation and stellar winds from th e cluster shape 



large gaseous pillars at the edge of the clou d ( Brand ner et al 



l2000h . CS observations bv lNurnberger etalj d2002l) show that 
the molecular gas also extends much furtherto the north (> 10') 
and to the south, hosting numerous massive molecular clumps 
which may be future sites of star formation. A number of 



prominent IR sources have been found so fa r dFrogel et al 



1977tlNurnberger)200 8: Nurnb ergeret alj2010l) as well as pro- 
plyd like objects dMiicke et alJl2002h . Rece nt AKARI observa- 



tions of NGC 3603 have been presented by Okad a et al. I d2010h 
showing that the [CII] 158 /mi emission is widely distributed 
and that the [OIII] 88 /mi emission follows the MIR, FIR, and 
radio continuum emission. 
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3. Observations 

We used the NANTEN2-4m antenna at 4865m altitude 
in Pampa la Bola in northern Chile to map the central 
2'x 2'region of NGC 3603. Our reference position is 
R.A.(J2000.0)=1 l h 15 m 08 s .85, Dec.(J2000.0)=-61°16'50'.'0, 
73" south-east of the central OB cluster NGC 3603 YC at 
R.A.(J2000.0)=ll h 15 m 07 s .26, Dec.(J2000.0)=-61°15'37/48. 
We observed the rotational transitions of 12 CO J=4-3 
(461.0408 GHz) and J=7-6 (806.6517 GHz) and the two 
fine-structure transitions of atomic carbon [Cl], 3 Pi- 3 Po 
(492.1607 GHz) and 3 P 2 - 3 Pi (809.3446 GHz) (henceforth 1-0 
and 2-1), between September and November 2006 with a dual- 
channel 460/8 10 GHz receiver. The exact line frequencies were 
taken fro m the Cologne Database f or Molecular Spectroscopy 



CDMS (MulleretalJ 12005, 2001). Double sideband (DSB) 



receiver temperatures were « 250 K in the lower frequency 
channel and « 750 K in the upper channel. The intermediate 
frequencies (IF) are 4 GHz and 1 .5 GHz, respectively. The lat- 
ter IF allows simultaneous observations of the CO 7-6 line in 
the lower sideband and the [C I] 2-1 line in the upper sideband. 
These two lines were observed simultaneously with one of the 
two lower lines in the 460 GHz channel. As backends, we used 
two acousto optical spectrometers (AOS) with bandwidths of 
1 GHz. The channel spacing was 0.37 km s _1 at 460 GHz and 
0.21 km s _1 at 806 GHz. The pointing accuracy was checked 
regularly on Jupiter, IRC+10216, and IRc2 in OrionA. The 
applied corrections were always <20" and usually < 10". To 
determine the atmospheric transmission, we measured the 
atmospheric emission at the reference position. Spectra of the 
two frequency bands were calibrated separately, and sideband 
imbalances were corrected using the atmospheric model atm 



(Atm ospheric Transmission at Microwaves, cf. iPardo et al 



20011) . Observations were taken on-the-fly (OTF), scanning in 
right ascension at a speed of 2.5"/sec and sampling every 10". 
The reference position was observed at the beginning of each 
OTF scanning line. 

The half power beam widths (HPBW) deconvolved from 
the observed full widths at half maximum (FWHM) are 
38.0"and 26.5"in the lower and upper receiver bands, respec- 
tively. Beam efficiencies fl e ff are 50% a nd 45%, respectively 



20081) . The median RMS 
and 1.25 K for CO 4-3, 



( Sim on et al.l 12007b iRramer et al 
across the map is 1.74, 1.53, 0.70 
7-6, [C I] 1-0, and 2-1 respectively. The forward efficiency, 
^eff = 86% in both bands, was determined from sky-dips. 
The raw data were calibrated to antenna temperatures T* A and 
scaled to main beam temperatures T m b with the factor F^/B^. 
We present all data in units of T,,,/,. Fifth order polynomial 
baselines were subtracted from all spectra. Observations of the 
atomic carbon lines are very challenging and require many OTF 
coverages to achieve a signal to noise (S/N) ratio above 3 <x. 



4. Data 

4.1. Dust and PAHs 

In Fig 0] we show a map of integrated CO 4-3 intensities from 
NANTEN2 overlayed to an 8 //m image taken with the Infrared 



Array Camera (IRAC) on board the Spitzer space telescopeQ. 
The position of the OB cluster is marked by a star. The im- 
age shows prominent molecular clum ps southwest and south - 
east of the cluster position. Following N iirnbereer et alj (2002), 
we refer to them as MM1 and MM2, respectively. The IRAC 
image shows the sharp PDR interfaces between the HII re- 
gion and the molecular cloud. The positions of these interfaces 
matc h the pillar-like struct ures visible in HST and VLT im- 
ages (Bra ndner et al. I l2000h . Note that these pillars are much 
smaller than the molecular clouds from which they protrude. 
We also note a dusty filament, which connects MM1 and MM2 
in a large arc 1.5' to 2' south of the western pillar head. This 
filament is not visible in CO or [CI] emission. The strong 
8 pm emission (> 10 4 MJy/sr) includes strong PAH emission 
(7.6,7.8, and 8.6 pm). The 8 pm band is an excellent tracer of 
the PDR interface. For a detaile d analysis of the emissio n by 
PAH's and very small grains see lLebouteiller et aL ( 20071) . 



4.2. NANTEN2 observations 
4.2.1 . Maps of integrated intensity 

Fig. [2] shows maps of integrated CO 4-3, 7-6, [C I] 1-0, and 
2-1 emission integrated over a velocity range from 4 to 
26 km s _1 . The two molecular clumps MM1 and MM2 are 
visible: MM1 to the southwest of the OB cluster and MM2 
to the southeast. The peak intensities are 267 Kkm s _1 for 
CO 4-3, 121 K-km s 1 for CO 7-6, 28 K-km s _1 for [C I] 1-0, 
and 30 Kkms" 1 for[Cl] 2-1. 

MM2 is very prominent in all 4 observed transitions while 
MM1 is weaker in the fine structure transitions. 

Both CO maps mark the transition from the HII region 
to the mol ecular cloud. MM1 is p resumably smaller than our 
beam size dNiirnberger et alj|2002l) . hence it is not possible to 
infer detailed structure information from the maps. 

The observed intensity distributions of the [C I] fine struc- 
ture transitions are more extended with a lower contrast in com- 
parison to the intensity distributions of the CO rotational line 
transitions. The peak positions of [CI] are slightly shifted away 
from the cluster compared to CO. However, the S/N of both 
[C I] transitions is too low to allow further conclusions. MM1 
can be identified in both CO maps. In [C I] a diffuse emission 
is visible without a clear distinction between clump center and 
its more diffuse environment. The MM2 peak positions of both 
CO maps match within the pointing accuracy. The same is true 
for both [C I] maps. Both [C I] peaks are shifted away from the 
OB cluster with respect to CO. However, we note a good spatial 
correlation between the CO and [C I] emission as all peaks co- 
incide within a beam radius. We do not see any C-CO layering 
as expected from simple edge-on, plane-parallel PDR scenar- 
ios. It is possible to explain coincident [C I] and CO emission 
with a face-on configuration, but the visible pillar structures in 
the HST observations and the very sharp interface visible in the 



1 Post-BCD data were retrieved from the Spitzer archive (URL 
http://archive.spitzer.caltech.edu/!. Because of strong saturation ef- 
fects, we chose the short exposure data. 
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Table 1. Observed integrated intensities I - J T m \, dv in units of K-km s 1 , Vlsr in units of km s and the line widths FWHM 
in units of km s computed as th , 1 st , and 2 nd moment of the spectra between v= 8. ..22 km s . The positions lie along a line 
starting from the OB cluster. The cluster position has been omitted. Rotational transitions of CO are denoted 43 ( I2 CO 4-3) and 
76 ( 12 CO 7-6). The atomic carbon fine-structure transitions are denoted 10 ( 3 P\ - 3 Pq) and 21 ( 3 f2 - 3 ^l)- 
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Fig. 2. Velocity integrated maps (200" x 200") of CO 4-3, 
7-6, [C I] 1-0, and 2-1 smoothed to a common angular res- 
olution of 38" (1.3-1.5 pc at 7-8 kpc), and integrated over 
a velocity range from 4 to 26 km s . Color scale and con- 
tours show the same data. Contours range between 10 and 90% 
of the peak intensities which are 267 K-km s for CO 4-3, 
121 K-km s _1 for CO 7-6, 28 K-km s" 1 for [Cl] 1-0, and 
30 K-km s _1 for [Cl] 2—1. The position of the compact OB 
cluster is marked by a white star. The position of the the pil- 
lars, as seen in the HST images, are marked by white squares. 
Dashed lines and small crosses spaced by 22" mark two cut 
from the OB cluster to the peak positions of the two clumps 
MM1 (short cut) and MM2 (long cut). The circle marks the 
resolution and t he two white crosses ma rk the position of the 
CS(2-1) peaks jNurnberger et alj l2002h in MM1 and MM2. 
The position (0,0) corresponds to R.A.(J2000.0)=1 l h 15 m 08 s .85, 
Dec.(J2000.0)=-61°16'50'.'0. 



IRAC image in FigQ]clearly show a prominent, edge-on PDR 
interface. 



The velocity structure of the observed field is shown in 
the CO 4-3 velocity channel maps (Fig. [3]). From the channel 
maps, it is evident that the emission towards MM1 shows a 
smaller kinematic range than in MM2. MM1 is not observable 
at velocities higher than 16 km s _1 , while the emission in MM2 
extends up to 21 km s _I . At larger velocities, the emission grad- 
ually shifts to the southwest. A similar kinematic behavior is 



also v isible from CS (3-2) observations by Nurnberger et al 



(2002). Their results also show a velocity drift for MM1 with 
distance from the OB cluster, possibly indicating a placement 
of MM1 in front of the cluster. In the appendix we show the 
velocity channel maps of the remaining 3 transitions. 




Fig. 3. Velocity structure of the observed field: CO 4-3 in- 
tensity in velocity channels of 1 km s _! width. Contours 
range between 5 and 45 K in steps of 5 K. The 
position (0,0) corresponds to R.A.(J2000.0)=1 l h 15 m 08^85, 
Dec.(J2000.0)=-61 o 16'5070. 



4.2.2. Spectra along a cut through the region 

To study the transition from the ionized to the molecular gas, 
we selected two cuts from the OB cluster to the peak positions 
of MM1 and MM2. The positions along the cut are 22" apart 
and marked in Fig. [2] The spectra along the two cuts, shown in 
Fig HI trace the kinematic structure of the interface regions and 
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Fig. 4. Spectra along two cuts connecting the OB cluster and 
the centre of MM1 (bottom) and MM2 (top). The CI spectra 
appear to show sinusoidal standing waves, however a detailed 
analysis shows no consistent standing wave pattern. These 
baseline variations are thus attributed to instabilities of the sys- 
tem and/or the atmosphere during the measurements. All data 
are at a common resolution of 38" and on the main beam tem- 
perature scale. 



the molecular clumps. The spectra at the positions 30"/- 10", 
40"/ -30", and 50"/-50 " do no t show simple Gaussian line pro- 
files. INiirnberger et al. I d2002l) observed similar line profiles in 
their CS 3-2 and 2-1 observations of MM2. They used two 
Gaussian components separated by 3 km s _1 to fit the emission 
profiles. Our analysis showed that it is not possible to use a sim- 
ilar two component Gaussian fit to reproduce our observations 
at all positions. We derived the first three moments of the spec- 
tral line, i.e., integrated intensity, mean velocity, and FWHM, 

6-22 km s _1 . The de- 



over the channels in the range v 
lived moments are given in Table Q] At the offsets -10"/70"and 
0'750"we don't find sufficient emission to derive significant 
values. 



4.3. Analysis 

4.3.1 . LTE - Column densities and temperatures 

In Table [2] we give the velocity integrated line ratios R74 
(hbl hi), R21 ihilho) and R i4 (/10//43) for all positions along 
the two cuts. We compute the LTE temperatures using the opti- 
cally thick CO 4-3 emission assuming a beam filling of 1 and 
obtain temperatures of 42 K for the peak position in MM2 and 
24 K for MM1. This is roughly consistent with estimates from 
the [C I] line ratio. R21 is a sensitive function of the [C I] ex- 
citation temperature. In the optically thin limit and assuming 
LTE, r ex = 38.3 K/ln[2.1 1/^21]. In MM2 we find ratios be- 
tween 0.8 and 1.1 which correspond to temperatures between 
40 and 60 K. For MM1 the [C I] line ratio gives a temperature 
estimate of 42 K, higher than the estimate derived from the CO 
4-3 emission. The difference can be explained by beam dilu- 
tion effects since MM1 is not resolved within the NANTEN2 
beam. However, the atomic carbon emission is heavily affected 
by noise. Another possibility is to use the LTE ratio of the two 



optically thick lines ^74 to derive an excitation temperature: 
T ex = 99 K/ ln[3. 06/^74]. Rj4 is remarkably constant across 
both clouds. We find an excitation temperature slightly above 
50 K for both clouds. First, this confirms that both clouds have 
comparable excitation conditions despite their different appear- 
ance. Second, they seem not to be composed of hot, unresolved 
clumps which would result in a much higher excitation tem- 
perature derived from Instead the gas appears to be more 
smoothly distributed and excited. 

The atomic carbon column density was derived under the 
assumption of LTE and optically thin [C I] emission. To de- 
rive the total H2 column densities we use the integr a ted in - 
tensity of C ls O 2-1 as presented in INiirnberger et alJ (120021) . 
Assuming opticall y thin emission, an isotope ratio of 500 and 
CO:H 2 =8xlO" 5 jFrerkinget al.ll 19821) . i.e., an assumed ratio of 
C 18 O/H2=2xl0~ 7 , we compute the column densities along the 
two cuts. The results are summarized in Table [2] For MM1, 
we find a total column density of 1.5 xlO 22 cirT 2 ; for MM2, 
N = 2... 5 x 10 22 crrT 2 . We find that the relative abundance of 
atomic carbon Nc/N tot drops towards the peak positions with 
values of 10~ 5 (MM1) and 8 - 19 x 10~ 6 (MM2), and in- 
creases at positions towards as well as away from the OB clus- 
ter. Even though the spatial resolution of our maps is rather 
coarse, this could be an indication for a C/CO stratification, be- 
cause in this case we would expect limb brightened [CI] emis- 
sion. The gas pressure within MM1 and MM2 (assuming a den- 
sity o f 10 4 cm" 3 ) is abou t 10-30% of the pressure in the HII re- 
gion ( Sha ver & Gossll970l) . Hence, the interface region is com- 
pressed and driven against the molecular cloud, visible, e.g., in 
the pillars protruding from MM1 and MM2. From the total col- 
umn density we can estimate the total mass per beam. The re- 
sults are given in Table|2] We also derived virial mass estimates 
across the cuts assuming a density of 10 4 cirT 3 and kinetic tem- 
peratures of 50 K (the densities are in agreement with the cal- 
culations from section l4~.3.3l l. The virial masses agree very well 
with the LTE masses and are much larger than the Jeans mass 
of the material of 87 M0. This means, the gas is gravitationally 
bound and not transient. We expect that the strong HII pres- 
sure and the gravitation al instability causes stro ng fragmenta- 
tion and cloud collapse. iNiirnberger et alJ d2002) derived a star 
formation efficiency of > 30% and a mean star formation rate 
of 1.3 x 1( T 3 M R /yr. Embed ded IR sources have already been 
identified (Fro g el et al 1 ll977h along the outskirts of the molec- 
ular clouds and particularl y at the base of the pillar like struc- 
tures at the edge of MM2 (Niirnberger 20081: INiirnberger et"al 
120101) . 



4.3.2. Dust Temperatures and FUV intensities 

To estimate the total FUV flux, we consider the total lumi- 
nosities of the most massive O stars in the central OB cluster. 



Recently, Mel ena et alJ (120081) published an updated census of 
the massive star content of NGC 3603. The most massive mem- 
bers of the OB cluster are 3 WNL stars, 14 03( and 3.5) stars 
(III+V) and at least 20 late O-type stars. They produce the FUV 
flux that dominates the heating of the PDRs and the molecular 
clouds. At an effective temperature of 24000 K, the ratio of 
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Table 2. Line ratios R74 (/76//43X R21 (hi/ ho) and R14 (/10//43) for all positions along the two cuts. The LTE excitation tempera- 
ture T ex ,43 is derived from the optically thick I43 emission, T ex ,74 is derived from R74 The C column densities are derived assuming 
LTE a nd optically thin [C I] emission. Values of the integrated intensities of C 18 2-1 are given as presented in Niirnberger et alJ 
(2002) along the two cuts, the LTE excitation temperature is derived fro m Z43 . The LTE colu mn densities of C ls O, and the total 
gas column density assume an isotope ratio of 500 and CO:H2=8xl0~ 5 dFrerking et alJH982l) . i.e. a ratio of C 18 O/H2=2xl0~ 7 . 
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Fig. 5. Left: Values of the FUV field in units of the Draine field derived from the HIRES 60, and 100 fim fluxes as color map. 
White, dashed contours mark the Spitzer/IRAC 5.8/im flux smoothed to a resolution of 35". The contour values are 25, 50, 
and 75% of the peak value. Overlayed in thin contours is the NANTEN2 CO 4-3 map. The central OB cluster is marked by 
a star. The two pillar heads are denoted by the squares, White crosses mark the positions along the cuts into MM1 and MM2. 
Right: Spitzer/IRAC 5.8 fim observation (in MJy/sr) overlayed with FUV fluxes as derived from IRAS/HIRES flux ratios. The 
contours correspond to values of x = 2000, 3000, 4000, and 5000. IRAC pixel size is 1.2". The position (0,0) in both panels 
corresponds to R.A.(J2000.0)=1 l h 15 m 08 s .85, Dec.(J2000.0)=-61 o 16'5070. 



the FUV energy density to the total energy density emitted by 
the star <I>fuv /^tot ~ 0.7 is maximal if we assume pure black 
body emission. At higher temperature, relatively more energy 
is emitted in the EUV range (13.6eV</iv< 130eVjl 

Applying effective temperatures an d luminosities a s given 
bv lMartins et al. I d2005h for t he O stars, Panagia (119731) for the 
B stars and CrowtherT i 20071) for the WR stars, we calculated 



a total FUV luminosity of 1.1 x 10 43 erg s _1 from the contri- 
bution of each cluster member. This is a lower limit since the 
spectral classes for more than 10 cluster members remains un- 
known. Depending on the distance from the cluster center d, 



2 IPauldrach et al] i 1998T) showed in model calculations that for ef- 
fective temperatures below 45000 K, the ratio of FUV to EUV photon 
rates is considerably high er than the r atio derived for pure blackbody 
spectra (see also Brandner et al. 200oj). This is because EUV photons 
are absorbed in the ionization front. We neglect this effect. 



in units of pc, the FUV flux is x - 9.4 x 10 4 oT 2 in units of 
the Draine field. At a distance of 7-8 kpc, one parsec corre- 
sponds to angular distances of 26-29", i.e., somewhat smaller 
than the projected distance between the cluster center and the 
edge of MM2 as seen for example in Fig.[T](«i 46"). If we as- 
sume that the cluster and the molecular cloud are situated in the 
same plane, the FUV field at the peak position in MM2 drops 
to x ~ 5000 - 6600. This is of the same order as the flux at 
the MM2-peak derived from the IRAS data. If the clumps are 
displaced from the plane of the OB cluster, or if the distance to 
the complex is different, the derived ;f changes. Uncertainties 
in the in-plane displacement of + 4 pc and in the distance (7-8 
kpc) lead \a X = 2400 - 6600. 

We use IRAS data to derive dust temperatures and the FUV 
continuum in the observed region. Enhanced resolution images 
of ~ 1' resolution were created using the maximum correlation 
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method ( Aumann et al. 1990). We estimate dust temperatures 
from the ratio of HIRES 60 fim and 100 /im data, assuming 
a dust spectral index of emissivity of 1.5. We obtained high 
resolution (HIRES) 60 and 100 fim images (2°x2°) from the 
IPAC data centei0. Along the two cuts, we find almost constant 
dust temperatures of 35 K. This value is very cl ose to the dust 
temp erature of 37 K derived from MSX data bv lWang & Chenl 
(2010). Again, comparison with the higher resolution Spitzer 
maps shows that the warm dust is not distributed homoge- 
neously. We conclude that the HIRES data are, therefore, af- 
fected by beam filling 
I - 



Following iNakapawa et al.l 4.1998), we combined 60 and 
100/im data to create a map of far-infrare d intensities between 
42.5 /im and 122.5 jm\ (IHelou et al 1 ll988l) . Under the assump- 
tion that all FUV energy absorbed by the grains is re-radiated 
in the far-infrared (FIR), we com pute the FUV fluxes (x) from 
the emergent FIR intensities (cf. iRramer et al.ll2005l) / fir, us 



ing y /yp = An 7 F i R , with^o = 2.7 x 10 3 erg s 1 cm 2 dDraine 



19781) . The resulting map of the FUV continuum is shown 
in the left panel of Fig. [5] Overlayed as thin contours is the 
NANTEN2 CO 4-3 map. At the peak position of MM1, we 
find FUV fields of x = 3460. Along the cut in MM2, we find 
that^f gradually increases from 3350 at the position of the OB 
cluster to 5230 at the peak position. Due to the low spatial res- 
olution of the HIRES data, these values are lower limits since 
beam dilution effects are not negligible. 

One would expect a decrease of FUV with growing dis- 
tance from the cluster. However, such a decrease is not seen, 
Fig. [5] shows the opposite. This is because the method of de- 
riving FUV intensities from IRAS data is not working for too 
small dust column densities because then UV radiation can not 
be effectively transformed to IR. Consequently, it would be best 
to exclude these parts from the further analysis. In our case, this 
is impossible due to the large beam. HIRES beam sizes vary 
from ~35"to over 120"across the maps, hence the FUV values 
in Fig. [5] are considerably beam diluted. In the right panel in 
Fig. [5] we overlayed the FUV fields derived from the HIRES 
data on 5.8fim Spitzer/IRAC observations of NGC 3603 to il- 
lustrate the beam filling effect: 60 and 100/mi emission can 
only be measured if dust is present at all. The Spitzer map 
shows no dust emission between the interface and the clus- 
ter. The HIRES beam size at the cluster position is almost cir- 
cular with a FWHM of ~ 35". Hence, a contribution to the 
IRAS/HIRES fluxes at the cluster position can only come from 
the interface regions of MM1 and MM2. If 60/100^m emission 
could be measured at l-2"resolution, one could not determine 
the FUV field towards the cluster and would observe a sharp 
interface between the H II region and the clumps. Instead we 
see no zero FUV intensity at the cluster and a gradual increase 
up to the peak 80" behind the interface. This is further illus- 
trated in the left panel in Fig.[5]where we additionally overlayed 
Spitzer/IRAC 5.8/im emission smoothed to the best HIRES res- 
olution of 35". At this resolution, up to 50% of the total flux is 
smeared across the H II region. 

At the MM2-peak, the HIRES beam is almost fully filled 
and the FUV fluxes are consistent with the assumption that the 



molecular cloud and the cluster lies in almost the same plane of 
the sky. This also supports the assumption that the FUV irradi- 
ation of the OB cluster created the pillars seen in the HST/VLT 
images. 

The derived FUV fluxes can be used to estimate the ratio 
^cii/^fuv, i-e- the fraction of FUV energy that is re-emitted in 
the the [CII] 158 /mi cooling line. We retrieved the [CII] spec- 
trurrQ from the ISO data archival and compute an integrated 
intensity of 2 x 10~ 3 erg s _I cirT 2 sr~'. Most of the FUV en- 
ergy that irradiates the cloud is re-radiated in the continuum, 
but a small fraction is re-emitted in cooling lines, such as the 
[CII] 158/im line. Applying the FUV flux derived from the 
HIRES data, we estimate /cii//fuv ~ 0.2. .0.3%. This is con- 
sistent with similar results from regions of active star forma- 
tion, e.g., nuclear regions of galaxies, where approximately 0. 1 
to 1 % of the total FUV input is re-radiated in the [CII] line 
dStacev etal.lll99ll) . However, this remains a rough estimate 
because of the large beam sizes involved. 

Considering the similar FUV estimates from stellar lumi- 
nosities and from the FIR data, we conclude that the FUV flux 
illuminating the molecular clumps MM1 and MM2 is 3000- 
6000 in units of the Draine field. Approximately 0.2 to 0.3% of 
the total FUV energy is re-emitted in the [CII] 158/vm cooling 
line. 

4.3.3. Escape Probability - density and 
temperatures 

We performed escap e probability (EP) calculations 

(Stut zki & Winnewisserl 1 19851) to derive local gas densi- 



ties, column densities, and gas temperatures from the observed 
line ratios (R74,R2i,Ri4, and Ic^oi-ilhi)- The model assumes 
a finite, spherical cloud geometry. We fitted the column 
densities of 12 CO and C independently. From the absolute 
CO 4-3 intensities we derived a filling factor I Q b s / I mo d- The 
total H2 column densities and masses have been calculated 
using CO:H 2 =8 x 10 5 . The results are shown in Table [3] 

For MM2, the temperatures from the EP calculations 
agree within 25% with the LTE temperatures for the positions 
30'7-10", 40"/-30", and 50'7-50". Masses and column densi- 
ties are comparable with results from the LTE calculations. The 
filling factor is roughly 0.3-0.6 across MM2. MM1 shows sig- 
nificantly smaller values of 0.1 - 0.2. This is consistent with 
the morphology namely that the clouds are clumpy and MM1 
is not resolved by the NANTEN2 beam. The derived temper- 
atures for MM1 are about a factor 2 larger compared to the 
temperatures derived from CO 4-3 and are compatible with 
the temperatures derived from the observed ratio R74 (consis- 
tent with / < 1). The gas densities found for MM1 and MM2 
are about 10 4 cm" 3 . In contrast to the LTE results, we find a 
monotonously growing Nc/Nh, along the cut in MM2 with val- 
ues of 1-5 xl0~ 5 . The column densities, given in Table [3] re- 
sult in C/CO abundance ratios between 0. 16-0.6 Similar abun- 
dance ratios have been observed in other Galactic star form- 
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Table 3. Results from the escape probability calculations. The filling factor / = 7 b s / I m od is derived from the CO 4-3 intensities. 
The CO and C column densities are results from the escape probability fit. 
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in g regions like Cepheus B, Orion or NGC7023 (; 
in IMookeriea et alil2006h . The C/CO ratio is su 



(see references 

IMookeriea et al.l 120061) . The C/CO ratio is supposed to de- 
pend on the local FUV intensity due to dependency of the 
CO formation and de struction balance on the FUV illumina- 
tion. I Sun et al.l d2008l) found significantly lower C/CO ratios in 



IC348 where the FUV field varies between 1 and 100 Draine 



units. On the other hand. lKramer et al.l (120081) found C/CO ra- 
tios comparable to our results in the Carina region, where the 
FUV field is also comparable to NGC 3603. We also find a 
monotonously growing C/CO abundance rat io along the cut 



in MM 2. A similar trend has been found by iMookeriea et al 



d2006h . namely that the C/CO ratio is higher further away from 
the sources of FUV radiation. 

5. Summary and Conclusions 

We used the NANTEN2-4m telescope to map the emis- 
sion of atomic carbon and CO in the vicinity of the central 
OB cluster in the Galactic star forming region NGC 3603. 
These data are the first observations of CO 4-3, CO 7-6, 
[CI] 1-0, and [CI] 2-1 in NGC 3603. We present fully sam- 
pled 200"x200"maps integrated over the full velocity range 
as well as velocity channel maps. The observed field includes 
the central OB cluster, as well as the two adjacent molecular 
clumps MM1 and MM2, hosting two pillars observed by the 
HST and similar in appearance to the famous pillar structures 
in the M16 nebula. We selected two cuts from the OB cluster 
position to the peak emission in MM1 and MM2 for detailed 
analysis. The spectra along the two cuts show a rich kinematic 
structure, especially towards MM2. The observed maps show 
a strong correlation between the spatial distribution of C and 
CO. This implies either a face-on configuration of the clouds, 
i.e., both clouds being not in the same plane as the OB cluster 
and thus being illuminated face-on, or a more complex con- 
figuration, for instance a composition of many small clumps. 
The overall good C-CO correlation may indicate an unresolved, 
clumpy structure, but the increase of Nc towards a cloud edge 
could hint towards relatively well shielded inner parts. 

Using escape probability model calculations we derived 
temperatures, densities, and column densities for both clouds 
from the observed line ratios. The temperatures are 43 K and 
47 K for the peak positions of MM1 and MM2, respectively. 
We find gas densities of n = 10 3 — 10 4 cm -3 in MM1 and 
MM2. From the best fit CO column densities we derive total 
H2 column densities (filling corrected) of 0.9 x 10 22 cirT 2 and 



2.5 x 10 22 cm" 2 for MM1 and MM2 and a ratio of N c /Nh 2 of 
l-5xl0~ 5 . 

The cluster is strongly interacting with the ambient molec- 
ular cloud and governing its structure and physical conditions. 
The stability analysis shows the existence of gravitationally 
bound gas which should lead to star formation. Embedded 
IR sources have already been observed in the outskirts of the 
molecular cloud and support our conclusions. 

We used HIRES/IRAS far-infrared data to narrow down 
the value of the FUV field at the positions of the two molec- 
ular clouds. Consistent with estimates from spectral type ap- 
proximations and with [CII] observations by ISO we find x ~ 
3 - 6 x 10 3 in units of the Draine field. 

However, many issues remain unresolved. For example, the 
analysis results are not conclusive regarding the dumpiness 
of the gas. NGC 3603 will be observed within the Herschel 
Guaranteed Time Key project Warm and Dense ISM - WADI. 
These data will allow us to perform a much more detailed study 
of the local gas conditions. Cooling lines like [O I] and [C II] 
as well as high-J CO lines can be used to study the dumpiness 
and the energy balance of the gas. 

Acknowledgements. We made use of the NASA/IPAC/IRAS/HiRES 
data reduction facilities. Data reduction of the spectral line data was 
done with the gildas software package supported at IRAM (see 
|http: //www, iram. f r/IRAMFR/GILDAS) . 

This work is financially supported in part by a Grant-in- Aid for 
Scientific Research from the Ministry of Education, Culture, Sports, 
Science and Technology of Japan (No. 15071203) and from JSPS (No. 
14102003 and No. 18684003), and by the JSPS core-to-core program 
(No. 17004). This work is also financially supported in part by the 
German Deutsche Forschungsgemeinschaft, DFG grants SFB494 and 
Os 177/1-1. 

References 

Aumann, H., Fowler, J., & Melnyk, M. 1990, AJ, 99, 1674 
Brandner, W., Grebel, E. K., Chu, Y.-H., et al. 2000, AJ, 119, 
292 

Crowther, P. A. 2007, ARA&A, 45, 177 
Crowther, P. A., Schnurr, O., Hirschi, R., et al. 2010, ArXiv 
e-prints 

Cubick, M., Stutzki, J., Ossenkopf, V, Kramer, C, & Rollig, 

M. 2008, A&A, 488, 623 
Dale, J. E.and Bonnell, I. A., Clarke, C. J., & Bate, M. 2005, 

MNRAS, 358,291 



Rollig et al.: Photon Dominated Regions in NGC 3603 



9 



Dedes, C, Rollig, M., Mookerjea, B., et al. 2010, A&A, sub- 
mitted 

Draine, B. T. 1978, ApJS, 36, 595 

Frerking, M. A., Langer, W. D., & Wilson, R. W. 1982, ApJ, 
262, 590 

Frogel, J. A., Persson, S. E., & Aaronson, M. 1977, ApJ, 213, 
723 

Gorti, U. & Hollenbach, D. 2002, ApJ, 573, 215 
Goss, W. M. & Radhakrishnan, V. 1969, Astrophys. Lett., 4, 
199 

Grabelsky, D. A., Cohen, R. S., Bronfman, L., & Thaddeus, P. 

1988, ApJ, 331, 181 
Harayama, Y., Eisenhauer, F., & Martins, F. 2008, ApJ, 675, 

1319 

Helou, G., Khan, I., Malek, L., & Boehmer, L. 1988, ApJ, 68, 
151 

Hollenbach, D. J. & Tielens, A. G. G. M. 1999, Reviews of 

Modern Physics, 71, 173 
Ikeda, M., Oka, T., Tatematsu, K., Sekimoto, Y., & Yamamoto, 

S. 2002, ApJS, 139,467 
Kramer, C, Cubick, M., Rollig, M., et al. 2008, A&A, 477, 547 
Kramer, C., Mookerjea, B., Bayet, E., et al. 2005, A&A, 441, 

961 

Lebouteiller, V., Bernard-Salas, J., Brandl, B., et al. 2008, ApJ, 
680, 398 

Lebouteiller, V., Brandl, B., Bernard-Salas, J., Devost, D., & 

Houck, J. R. 2007, ApJ, 665, 390 
Martins, F, Schaerer, D., & Hillier, D. J. 2005, A&A, 436, 

1049 

Melena, N. W., Massey, P., Morrell, N. I., & Zangari, A. M. 

2008, ApJ, 135,878 
Mookerjea, B., Kramer, C, Rollig, M., & Masur, M. 2006, 

A&A, 456, 235 

Miicke, A., Koribalski, B. S., Moffat, A. F. J., Corcoran, M. F, 

& Stevens, I. R. 2002, ApJ, 571, 366 
Muller, H. S. P., Schloder, F, Stutzki, J., & Winnewisser, G. 

2005, Journal of Molecular Structure, 742, 215 
Muller, H. S. P., Thorwirth, S., Roth, D. A., & Winnewisser, G. 

2001, A&A, 370, L49 
Nakagawa, T., Yui, Y. Y, & Doi, Y, e. a. 1998, ApJS, 1 15, 259 
Niirnberger, D. E. A. 2008, Journal of Physics: Conference 

Series, 131 

Niirnberger, D. E. A., Bronfman, L., Yorke, H. W., & 

Zinnecker, H. 2002, A&A, 394, 253 
Niirnberger, D. E. A., Vehoff, S., Hummel, C. A., & Duschl, 

W. J. 2010, Rev. Mex. Astron. Astrofis., 38, 71 
Okada, Y, Kawada, M., Murakami, N., et al. 2010, A&A, 514, 

A13+ 

Ossenkopf, V., Rollig, M., Simon, R., et al. 2010, ArXiv e- 
prints 

Panagia, N. 1973, AJ, 78, 929 

Pandey, A. K, Ogura, K, & Sekiguchi, K. 2000, PASJ, 52, 847 
Pardo, J., Cernicharo, J., & Serabyn, E. 2001, IEEE 

Transactions on Antennas and Propagation, 49, 1683 
Pauldrach, A. W. A., Lennon, M., Hoffmann, T. L., et al. 1998, 

in Astronomical Society of the Pacific Conference Series, 

Vol. 131, Properties of Hot Luminous Stars, ed. I. Howarth, 

258-+ 



Rochau, B., Brandner, W., Stolte, A., et al. 2010, ApJ, 716, L90 
Rollig, M., Abel, N. P., Bell, T., et al. 2007, A&A, 467, 187 
Rollig, M., Ossenkopf, V., Jeyakumar, S., Stutzki, J., & 

Sternberg, A. 2006, A&A, 451, 917 
Sakai, J. I., Nagasugi, Y, Saito, S., & Kaufmann, P. 2006, 

A&A, 457,313 

Shaver, P. A. & Goss, W. M. 1970, Australian Journal of 

Physics Astrophysical Supplement, 14, 133 
Simon, R., Graf, U., Kramer, C, Stutzki, J., & Onishi, T. 2007, 

NANTEN technical report, 13.2.2007, 1, 1 
Stacey, G. J., Geis, N., Genzel, R., et al. 1991, ApJ, 373, 423 
Stolte, A., Brandner, W., Brandl, B., & Zinnecker, H. 2006, AJ, 

132, 253 

Stolte, A., Brandner, W., Brandl, B., Zinnecker, H., & Grebel, 

E. K. 2004, AJ, 128, 765 
Storzer, H., Stutzki, J., & Sternberg, A. 1996, A&A, 310, 592 
Stutzki, J. & Winnewisser, G. 1985, A&A, 144, 13 
Sun, K, Ossenkopf, V., Kramer, C, et al. 2008, A&A, 489, 207 
Sung, H. & Bessell, M. S. 2004, AJ, 127, 1014 
Tatematsu, K, Jaffe, D. T, Plume, R., Evans, II, N. J., & Keene, 

J. 1999, ApJ, 526, 295 
Wang, J. & Chen, Y. 2010, Science in China G: Physics and 

Astronomy, 53, 271 
Zinnecker, H. & Yorke, H. W. 2007, ARA&A, 45, 481 



Appendix A: Velocity Channel maps 

Here we present the velocity channel maps of all observed tran- 
sitions. Each map is shown at its original spatial resolution, i.e. 
38" for CO 4-3 and C I 1-0 as well as 26.5" for CO 7-6 and 
CI 2-1. 
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Fig.A.l. Velocity structure of the observed field: CO 4-3 intensity in velocity channels of 1 km s _1 width. Contours range 
between 5 and 45 K in steps of 5 K. The position (0,0) corresponds to R.A.(J2000.0)=1 l h 15 m 08 s :85, Dec.(J2000.0)=-6ri6'5070. 
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Fig.A.2. Velocity structure of the observed field: CO 7-6 intensity in velocity channels of 1 km s _1 width. Contours range 
between 2 and 20 K in steps of 2 K. The position (0,0) corresponds to R.A.(J2000.0)=1 l h 15 m 08 s :85, Dec.(J2000.0)=-6ri6'5070. 
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Fig.A.3. Velocity structure of the observed field: [CI] 1-0 intensity in velocity channels of 1 km s _1 width. Contours range 
between 1 and 10 K in steps of 1 K. The position (0,0) corresponds to R.A.(J2000.0)=1 l h 15 m 08 s :85, Dec.(J2000.0)=-6ri6'5070. 




Fig.A.4. Velocity structure of the observed field: [CI] 2-1 intensity in velocity channels of 1 km s _1 width. Contours range 
between 1 and 10 K in steps of 1 K. The position (0,0) corresponds to R.A.(J2000.0)=1 l h 15 m 08 s :85, Dec.(J2000.0)=-6ri6'5070. 



